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Resonant behavior of the Raman process is predicted when the laser frequency is close to the 
orbital excitation energy of LaMnOa at 2 eV. The incident photon creates a vibrationally excited 
self-trapped "orbiton" state from the orbitally-ordered Jahn- Teller (JT) ground state. Trapping 
occurs by local oxygen rearrangement. Then the Franck-Condon mechanism activates multiphonon 
Raman scattering. The amplitude of the n-phonon process is first order in the electron-phonon 
coupling g. The resonance occurs via a dipole forbidden d to d transition. We previously suggested 
that this transition (also seen in optical reflectivity) becomes allowed because of asymmetric oxygen 
fluctuations. Here we calculate the magnitude of the corresponding matrix element using local spin- 
density functional theory. This calculation agrees to better than a factor of two with our previous 
value extracted from experiment. This allows us to calculate the absolute value of the Raman 
tensor for multiphonon scattering. Observation of this effect would be a direct confirmation of the 
importance of the JT electron-phonon term and the presence of self-trapped orbital excitons, or 
"orbitons." 



I. INTRODUCTION 

Manganese oxide materials attract attention because of the "colossal magnetoresistance" (CMR) phenomenon [0, 
and because of a very rich phase diagram [^| of ground states with competing order parameters. The Mn+^ ion of 
the parent LaMnOa compound has (^ig'^g) configuration with an inert t2g core (spin 3/2). The half-filled doubly 
degenerate eg orbitals {d^2_y2, d^^-2_j.2) are Jahn- Teller (JT) unstable. A symmetry-breaking oxygen distortion 
(resulting in the Mn-0 bond lengths of 1.91, 2.18 and 1.97 A lowers the energy of the occupied orbital. The 
corresponding orbitally ordered state sets in at Tjt = 750 K with x— and y— oriented Cg orbitals in the x — y plane 
with wavevector Q = (7r,7r, 0) The orbital order drives antiferromagnetic (A- type) spin order below the Neel 
temperature Tjy — 140 K. 

There is still controversy about the origin of the orbitally ordered state. Strong electron-electron correlations may 
lead to orbital order via the superexchange interaction Q which lifts degeneracy of the eg states 0. In another 
scenario proposed by Millis Q the Jahn- Teller electron-phonon (e-ph) interaction g j|] causes the orbital order and 
contributes to CMR. Extensive numerical work by Dagotto et al. showed that the two approaches give qualitatively 
similar answers. We prefer a model where the JT interaction g plus large Hubbard U and Hund energy Jh leads to 
single occupancy of the Mn eg levels and a gap to on-site d to d excitations, rather than assigning the gap purely to 
Coulomb interactions as in a multi-orbital Hubbard model. We believe that the importance of the JT interaction is 
evident at low hole doping of Lai_2;Sr(Ca)a;Mn03, whose insulating nature is naturally explained by formation of the 
anti- Jahn- Teller polarons [pd| . 

In this paper we present a detailed prediction of resonant multiphonon Raman features, whose observation would 
be a direct measure of the importance of the JT electron-phonon term. When the oscillator potential curves of 
ground and excited states are displaced relative to each other, then vibrational Raman scattering is activated by a 
; ^ Franck-Condon (FC) two-step mechanism. Our Hamiltonian for LaMnOa, with [/ — > cxd, leads to a picture where the 
, ground state and low-lying excited states are simple products of localized orbitals, one per atom. In the first step of 
the FC Raman process, the incident photon creates an orbital defect in the ordered JT ground state (one Mn ion has 
the upper rather than the lower state of the JT doublet occupied.) This Frenkel exciton (also called an "orbiton") is 
self-trapped |l2[ by oxygen rearrangement from the JT state. The FC principle has the oxygen positions undistorted 
during optical excitation, producing a vibrationally excited state of the orbiton. In the second step of the Raman 
process, this virtual excitation decays back to the orbital ground state, but not necessarily the vibrational ground 
state. The amplitude for ending in a vibrationally excited state is determined by displaced-oscillator overlap integrals. 
This allows n-phonon resonant Raman scattering with amplitude proportional to the first order of the e-ph interaction 
p3[ . The process is illustrated in Fig. |l]. In the conventional Raman scattering process, where electronically excited 
states do not alter atomic positions, the amplitude of the n-phonon peak is proportional to the nth order of the e-ph 
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interaction which is smaller by n — 1 orders of magnitude. The conventional process can be divided into three steps. 
(1) The incident photon creates an electron-hole pair (or exciton). (2) This electron-hole pair is scattered into another 
state by sequential emission of n phonons via n powers of the e-ph interaction 7ic-ph oc g. Higher-order interactions, 
such as the electron- (two phonon) interaction also enter, but do not increase the order of magnitude of the process. (3) 
The electron-hole pair recombines, emitting a scattered photon. In this formulation, the intensity of the two-phonon 
Raman process is smaller than one-phonon by several orders of magnitude 10~^-10~^, determined by the 2n-th power 
of the ratio of c-ph to electronic energies. 



II. FRANCK-CONDON MECHANISM 



We use a model Hamiltonian |11|, essentially the same as used by Millis [Q, with two eg orbitals per Mn atom, 
fully respecting the symmetries of the orbitals and the crystal. The electron-phonon term TYjt stabilizes the orbitally- 
ordered ground state via a cooperative JT distortion. Oxygen displacements along Mn-O-Mn bonds are modeled by 
local Einstein oscillators: 



(1) 



The interaction TYjt consists of linear energy reduction of an occupied d2^2_^2 orbital (the corresponding creation 
operator is cj, and number operator is fig^x = cj,(£)ca;(£)) if the two oxygens in the ztx direction expand outwards. 
Similar terms are included for y and z oxygens if d^y2_j,2 or c?3^2„r2 orbitals are occupied. The strength g—l.SA eV/A 
of the JT coupling g determines the JT splitting of the orbitals 2 A = 1.9 eV, and was chosen to agree with the lowest 
optical conductivity peak p^ ]. The displacement Ui^a is measured from the cubic perovskite position of the nearest 
oxygen in the d-direction to the Mn atom at £. The oxygen vibrational energy Tioj — hy/K/M = 0.075 eV is taken 
from a Raman experiment ]l6[ |. In addition there is an on site Coulomb repulsion U and a large Hund energy Jh- 
In the limit U ^ oo and Jh oo, electronic motion at half- filling is suppressed due to single occupancy of the Mn 
sites; additional orbital splitting caused by superexchange interactions is left out. 
The Hamiltonian H = Hjt + Ti-L gives an orbitally ordered ground state: 



io,o>=n4wn^i'(^')iw 



>, 



(2) 



where {0} refers to the lattice vibrational ground state with oxygen atoms in distorted equilibrium positions (Van 
Vleck Q2-type distortions) U£±x — ±uo, ui±y — ^uq, Ui±z = if £ A sublattice {exjp{iQ ■£) = !) and opposite sign 
distortions if £ G B sublattice {exp{iQ -f) = —1.) The magnitude of the distortion 2uq = •\/2A/Mcj 2=0.296 A agrees 
with neutron diffraction data within 10-15%, confirming the internal consistency of the model. Operators y 
create electrons with orbitals ^ x,y — —{d^z^-r^ T d,j;2_y2)/y/2 alternating on A and B sublattices. The lowest- lying 
electronic excitation of the Hamiltonian (|l|) is a self-trapped exciton or orbiton , which gives a broad line in the 
optical conductivity |^ centered at 2A « 2 eV. 

In the ri-phonon Raman process, incident light of frequency ujl is scattered with a shifted frequency los — i^l — nui. 
The Raman cross-section tensor R" can be found as following : 



{m},z 



< 0, n\eL ■ p|i, m >< i, 'rn\es ■ p|0, > 
A + N{m\huj — hujL + ijm 



NRT 



S{ujji — nuj) 



(3) 



where (Jq = is a Compton cross section (r^ = e^/rrigC^). The summation goes over all electronic states i and all the 
corresponding vibrational quanta {m}. Nlm} = mx + my + mz + m-x + m-y + m-z is the total number of vibrational 
quanta. The nonresonant term (NRT) is obtained from the resonant term by permuting el with is and changing 
—ujl to Lus- The imaginary frequency 7m takes into account the finite lifetime of the intermediate state \i,m >. The 
final state |0,n > has an electronic ground state plus n vibrational quanta. Summation over all the possible states 
with a total number of n vibrations is assumed in Eq. (^. For example, one phonon can be excited on any of the 6 
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neighboring oxygen atoms; two phonons can be excited in 21 ways; three phonons in 56 ways, and so on. The ground 
state couples to excited electronic states by the electron-radiation Hamiltonian (p- A). In LaMnOs, we consider only 
the lowest excited electronic state with an orbital flip (e.g. \X >-type to \Y >-type By neglecting coupling to 

higher electronic states we underestimate (perhaps by a significant factor) the first order Raman peak intensity. For 
multiphonon Raman scattering, in first approximation, we assume that only orbiton intermediate states contribute. 




FIG. 1. Schematic Franck-Condon mechanism for the multiphonon Raman process via the orbiton intermediate state. The 
lowest energy configuration of the orbiton has energy A and large oxygen distortions from the JT ground state (GS). The most 
probable intermediate state (the strongest resonance of the Raman process) occurs at ojl ~ 2A « 2eV, rather than ul ~ A. 



III. LSDA DIPOLE MATRIX ELEMENTS 

To evaluate dipole matrix elements in Eq. (^) we use the FC approximation. The wavefunctions |0, >, \i,m > 
and |0, n > are written as products of vibrational wavefunctions x(^) dependent on the oxygen positions i?, and 
electronic wavefunctions ip[r^R) dependent on both electronic r and vibrational coordinates. The electronic dipole 
matrix element is a d to d transition and therefore forbidden when the surroundings are symmetric. Searching for a 
mechanism to activate this transition, wc notice that an asymmetric oxygen displacement will cause Mn Cg orbitals 
to acquire an admixture of 4p character. A typical mixing coefficient is 

f dV 

Jz= drip3z2_^2-—ipz/{ed-ep) (4) 
J ouz 

where "02 is an orbital of p character, and dV/duz is the perturbation caused by a displacement of oxygen £ + z in z 
direction. The corresponding allowed optical matrix element is 



dfll^Sz^^r^Pzll^z- (5) 



The resulting dipole matrix element is 



< i, m\esx ■ p\0, n >= V" ^adaea < m,xm-xmym-ymzm-z\{ua + U-a)\nxn-xnyn-ynzn-z > (6) 



a—x,y,z 



If the ground state is described by Eq. (^), then from symmetry one can show that ^xdx — lydy — —^zdzl'^- 
In our previous work p2| a phenomenological parameter 7c? = —^zdz was introduced to account for the observed 
spectral weight of the optical conductivity peak due to the self-trapped exciton. The oscillator strength / defined as 
/ duja{uj) — (ttNc^ /2mefl)f is equal in our model to fzz — 2{{^d)'^ /meMu;'^)/{2A/huj + 1), and fxx — fyy = fzz/4:- 
Here N/ft is the Mn atom concentration. The measured spectral weight 540r2~^cm~^ eV of the lowest broad line 
centered at 2 eV [|l^ corresponds to /oxp — 0.113 or "fd — 1.7. Here we use density functional theory (DFT) to 
calculate an induced dipole matrix elements to test whether our choice of the parameter jd from the optical data was 
justified. 
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LaMnOs has been extensively studied by first-principles approaches ||T^ including the local-spin-density approxi- 
mation (LSDA) of DFT, LDA-I-U, and Hartree-Fock methods. Information about electronic and magnetic structure, 
and about electron-phonon and Coulomb interactions has been obtained. Here we use LSDA to calculate the dipole 
matrix element for d to d transitions as it is induced by asymmetric oxygen distortions. Rather than calculating 
dV/duz and doing perturbation theory as in Eqs. (^,^), we directly calculate < i\e ■ p\0 > in the presence of an 
imposed asymmetric oxygen distortion. 

To solve the LSDA equations we use the plane- wave pseudopotential method [ p^|p^ with a spin-dependent exchange- 
correlation potential and a supercell approach. Calculations were done for a 10-atom perovskite supercell with 
only (32-type JT oxygen distortions (the rotation of the MnOg octahedra is omitted.) The point group symmetry of 
the cell is thus D^h- For convenience, the magnetic order is taken to be ferromagnetic with 4.00 /is spin magnetization 
per formula unit. The lattice constant 3.936 A gives the same cell volume as observed for LaMnOs. The magnitude 
of the in-plane oxygen distortions along the Mn-O-Mn bonds is uq — 0.14 A. The self-consistent charge density was 
calculated using six special k points in the irreducible wedge of the Brillouin zone. Then the self-consistent potential 
was used to calculate wavefunctions at the T point. 
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FIG. 2. The LSDA wavefunctions ^'(r) at fc = (0,0,0) point for the Mn eg states Big (occupied) and Aig (empty). Arrows 
indicate the direction of the Q2-type oxygen displacements in the xy plane. Negative (dashed) and positive (solid) contour 
values are in units of 0.1 e^^^ /{a.u.)^^'^ 



A symmetry analysis of the pseudowavefunctions for k = (0, 0, 0) and plots of ji/jp around the Mn site allow us to 
distinguish Mn Cg states from other states. The eg {X- and F-type) orbitals form states of Aig and Big symmetries. 
Two of the four Mn Cg states are shown on Fig. ||. By introducing a small displacement of all the apical oxygens in +z 
direction along the Mn-O-Mn bonds, one induces a,pz dipole matrix element {Aig to Aig and Big to Big). Similarly, 
in-plane oxygen displacement in the +x direction induces a px matrix element. In table || we present LSDA results 
for these matrices. The imposed displacement has lowered the D^h symmetry of the supercell, permitting transitions 
between states below and above the Fermi level which were previously labeled as Aig and Big. The energies of the 
Mn Cg states do not alter much (less than 0.04 eV) for small distortions (0.016 A) and the induced dipole matrix 
elements are linear with oxygen distortion. 

It is convenient to measure the induced dipoles in units of {rrieMuj^y^^ — 0.472?i/(aB)^, where me and M are 
electron and oxygen masses and ae is the Bohr radius. The calculated dipoles {jd)z — 2.39 and {'yd)x^y = 1.9 give 
the oscillator strengths fzz ~ 0.43 and fxx,{yy) — 0.27. These answers are somewhat larger than experiment, but 
well within the expected accuracy of our model. This accuracy is limited by four factors: (1) convergence of the 
LSDA result; (2) neglect of rotational distortions; (3) simplification of magnetic structure to ferromagnetic; and (4) 
applicability of LSDA to strongly-correlated electrons. We tested the first by varying the number of plane waves used 
in the pseudopotential expansion. The answers reported here used a plane-wave cut-off -Epw = 135 Ry. The value of 
^d increases by 9% at E-p^ — 120 Ry and decreases by 13% at E-p^ = 100 Ry, so the convergence error is estimated 
at 10%. Sources (2) and (3) will causes errors of similar size, we believe. The biggest uncertainty is source (4). In an 
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unrelated DFT study on the self-trapped-exciton in crystalline NaCl |2^, we found evidence that LSDA incorrectly 
shifts the energy of localized states relative to delocalized states. However, the wavefunctions under consideration 
here are quite well-localized and agree qualitatively with expectation. Therefore we do not believe that point (4) is 
the source of a major error. Therefore we believe that this calculation provides firm evidence that phonon-activation 
is strong enough to account for the observation of d to c? (orbiton) transitions in optical reflectivity, as we assumed in 
our previous work, and as we now use to predict Raman spectra. 

IV. THE RAMAN TENSOR 

The absolute cross sections for the Raman process Eq. (0) can be evaluated using expression (B|): 



"'^ = o-Q— f (7arfa7/3rf/3)^^(^i? " nuj)'^S{n - N{f}) 



dujftdQ Lu , 

{/} 



y ncoA^,{m{f}) ^^^^ 

m— 



(7) 



{/}) - ^ <5(m - N{m'}) < f\u^ + u 

— a. 

{m'} 

where the induced dipole matrix elements 'yd and displacements u are measured in units of (TOgMo;^)"'^/^ and yJh/Muj 
respectively. In order to evaluate vibrational overlap integrals AQ_/3(rn, {/}) one needs the expressions for overlap 
integrals of displaced harmonic oscillators. 



<niK>=(-ir""^y;^e--/^A£"^-"^^(-l)'=— — -— , i/ni>n2 (9) 

^ k\{n2 — ky\ni — n2 + k)\ 

, , K / no — ni\ , , . 

< ni\u\n2 >= — ( IH — ) < ni\n2 >, (10) 



^/2 



where k is related to the Jahn- Teller gap as A = Ak^Hlu. The overlap < ni|n2 > for ni < n2 has the same expression 
as with ni and n2 interchanged and the sign of the displacement uq = y/2K changed (k — > —k). When using 
expressions (^, ^0|), the signs of the i + x, i — y oxygens displacements are positive and i — x, £ + y are negative, if 
£ €: A (reverse signs ii £ G B) and £zL z oxygens are undisplaced. Evaluation of the overlap integrals Aap^m, {/}) is 
straightforward. For example, for the first order Raman peak, only four one-phonon final states will contribute: 

-A Am 

Aa/3(m, {/} = 1 . . .4) = <5„,^ + 1 - m) (11) 

ml 4kA 

with no contribution to the non-diagonal part of the tensor a ^ j3. For second- and third-order Raman scattering, 
four and eight final states contribute to the non-diagonal part of the tensor: 

Aa:v{vz,zx)[m,{f} = 1 . . A) = ± — , 

Axy(y.,.x){m, {/} = 1 . . . 8) - ±_____(A + l^Tn) (12) 

The formulas for the diagonal parts of the higher-order Raman tensors are more complicated and will not be given 
here. 

To model the damping term 7^ of vibrational level m, we use expression 7™ — 7o\An + li as in a sequence of 
convolved Gaussians, intended to mimic the local densities of phonon states on oxygen atoms. The value 70 = 120 
cm~^ was taken. The Raman cross section shown on Fig. ^ has a pronounced resonant behavior when the laser 
frequency ojl approaches the orbiton energy 2A. The first-order cross-section as seen in Fig. ^ is underestimated, 
because we neglected coupling to higher electronic levels in Eq. (^, whose contribution to multiphonon peaks is 
probably negligible. The polarization dependence of the cross-section is only due to the dipole matrix element effect, 
namely Tl]:^''^ = Tiyy'^ = 7?.^f '^/16, R'^'^yz = Rl'^ = ^^Rly- The anisotropy of the optical conductivity a^x/cTzz is 
quadratic and Raman intensity TZxx/T^zz fourth-power in the dipole matrix-element anisotropy "fxdx/jzdz- Actual 
occupied Cg orbitals may be rotations of our idealized state Eq. (^) in the Cg space. Even small deviations from Eq. 
(0) might cause a noticeable change of the ratio "fxdxllzdz from the value 0.5. Therefore the predicted anisotropics 
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of 4 and 16 for optical and Raman spectra are not necessarily robust, but the relation Hxx/Ti-zz = i^xxl'^zzf' should 
hold. 

Published Raman measurements | [l6t on undoped LaMnOs do not extend to the multiphonon region and resonant 
behavior has not been tested experimentally. Most experiments use the Ar"^ laser, for which frequency = 2.41 eV, 
we predict the multiphonon cross sections Tl^^ = 0.514, Tl\^ = 0.243, Tl\^ = 0.056, Tl\^ = 0.025 in units of cto sr'^ 
Recently some features in Raman spectra on LaMnOa around 1100 cm~^ were reported by several groups p3| , p^ . 
These are probably the effect we are predicting. 



2.0 




Laser Frequency 0)^ (eV) 

FIG. 3. The absolute value of the multiphonon Raman cross section per solid angle versus incident photon energy. Resonant 
behavior is predicted for ujl close to orbiton energy 2A = 1.9 eV. The damping constant is 70 — 120 cm~^, the induced dipoles 
are —jzdz = 2'yxdx = 1.7 (meMuj'^)^^^ , cross section unit is ao = (e^/meC'^)^. 



V. CONCLUSION 



We advocate a picture of the orbitally ordered state of LaMnOa where electron-phonon interactions (in a context of 
large Hund and Hubbard energies) have a major influence. Our picture is disputed by other theorists |2^]. Therefore 
we attempt here to provide predictions which can qualitatively distinguish our model from others. In our model, the 
lowest electronic excitation is the 1.9 eV transition across the Jahn- Teller gap, modified by self-trapping to give a 
minimum gap half as large. This has successfully described the observed optical gap as a Franck-Condon broadened 
self-trapped exciton. The present paper uses density-functional theory to eliminate the need for a phenomenological 
coupling 7d to account for this transition. 

As a more stringent test, we here predict a new feature unique to the FC physics of the self-trapped exciton, namely 
a sequence of resonant multiphonon Raman peaks. We predict the absolute values of the multiphonon Raman cross 
section tensors. 

A hot luminescence process can also give rise to a multiphonon peaks. Incident light can excite an orbiton with a 
long lifetime which can recombine after vibrational energy nuj being lost in the intermediate state through anharmonic 
interaction. The question whether the Raman or hot luminescence mechanisms dominate the scattering intensity has 
an old history . In the Raman case the intensity of the higher order peaks should decrease and the width of the 
peak representing a convolved density of local phonon modes should increase. In the hot luminescence the intensities 
of the higher order peaks are of the same order and the width of the lines decrease with increasing order. In addition 
a strong emission peak can be observed at the exciton absorption edge. Raman techniques can serve as a direct probe 
of the orbiton excitation in LaMnOa. 
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TABLE I. The absolute values of the induced dipole ma- 
trix elements pz (left) and px (right) per unit oxygen displace- 
ments in units of (meMo;^)^^^ between the LSDA wavofunc- 
tions of the Ajg and B^g symmetry before an additional asym- 
metric oxygen distortion. The resulting oscillator strengths 
for the transition between the occupied (occ) and empty (em) 
states are fzz = 0.43 and fxx,(yy) — 0-27, which corresponds 
to (7d)^ = 2.39 and i'yd)x,y = 1.9. 





A occ 

^Ig 


A em 


ryocc 

Big 


Tjem 

Big 


A OCC 

^Ig 


A em 
^Ig 


T-jOCC 

-Dig 


rjcm 
-Dig 


A OCC 
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3.37 











1.83 


2.03 


0.90 


A em 
^Ig 


3.37 











1.83 





1.55 


0.05 


-ryocc 

Big 











0.17 


2.03 


1.55 





0.82 


oem 








0.17 





0.90 


0.05 


0.82 
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